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Abstract. We examine the algebraic complete integrability of Lotka-Volterra 
equations in three dimensions. We restrict our attention to Lotka-Volterra 
systems defined by a skew symmetric matrix. We obtain a complete classifi- 
cation of such systems. The classification is obtained using Painleve analysis 
and more specifically by the use of Kowalevski exponents. The imposition 
of certain integrality conditions on the Kowalevski exponents gives necessary 
conditions for the algebraic integrability of the corresponding systems. We 
also show that the conditions are sufficient. 



The Lotka-Volterra model is a basic model of predator-prey interactions. The 
model was developed independently by Alfred Lotka (1925), and Vito Voltcrra 
(1926). It forms the basis for many models used today in the analysis of population 
dynamics. In three dimensions it describes the dynamics of a biological system 
where three species interact. 

The most general form of Lotka-Volterra equations is 



We consider Lotka-Volterra equations without linear terms (e^ — 0), and where 
the matrix of interaction coefficients A = (aij) is skew-symmetric. This is a natural 
assumption related to the principle that crowding inhibits growth. The special 
case of Kac-van Moerbeke system (KM-system) was used to describe population 
evolution in a hierarchical system of competing individuals. The KM-system has 
close connection with the Toda lattice. The Lotka-Volterra equations were studied 
by many authors in its various aspects, e.g. complete integrability [7 Poisson and 
bi-Hamiltonian formulation ([l^, [H], [E]), stability of solutions and Darboux 
polynomials (0, [T5j). 

In this paper we examine the algebraic complete integrability of such Lotka- 
Volterra equations in three dimensions. The basic tools for the required classi- 
fication are, the use of Painleve analysis, the examination of the eigenvalues of 
the Kowalevski matrix and other standard Lax pair and Poisson techniques. The 
Kowalevski exponents are useful in establishing integrability or non-integrability of 
Hamiltonian systems; see [T], [2], [5], [T^, [H], [T7], [U]. The first step is to impose 
certain conditions on the exponents, i.e., we require that all the Kowalevski expo- 
nents be integers. This gives a finite list of values of the parameters satisfying such 
conditions. This step requires some elementary number theoretic techniques as is 
usual with such type of classification. In the three-dimensional case the general 
expressions for the Kowalevski exponents are rational and therefore the number 
theoretic analysis is manageable. 



1. Introduction 
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The second step is to check that the leading behavior of the Laurent series 
solutions agrees with the weights of the corresponding homogeneous vector field 
defining the dynamical system. In our case the weights are all equal to one and 
therefore we must exclude the possibility that some of the Laurent series have 
leading terms with poles of order greater than one. This is a step usually omitted 
by some authors due to its complexity, but in this paper we analyse this in detail. 
To accomplish this step we use old-fashioned Painleve Analysis, i.e., Laurent series. 
The application of Painleve analysis and especially of the ARS algorithm (see [3] 
0, [2], [8], [14]) is useful in calculating the Laurent solution of a system and 
check if there are {n — 1) free parameters. 

In performing Painleve analysis we use the fact that the sum of the variables 
is always a first integral. Surprisingly the Painleve analysis does not reveal any 
additional cases besides the ones already found by using the Kowalevski exponents. 
In this classification of the algebraic completely integrable Lotka-Volterra systems 
we discover, as expected, some well known integrable systems like the open and 
periodic Kac-van Moerbeke systems. 

To make sure that our conditions are not only necessary but also sufficient we verify 
that the systems obtained are indeed algebraically completely integrable by checking 
the number of free parameters. We also have to point out that our classification is 
up to isomorphism. In other words, if one system is obtained from another by an 
invertible linear change of variables, we do not consider them as different. Modulo 
this identification we obtain only six classes of solutions. 

The Lotka-Volterra system can be expressed in hamiltonian form as follows: 
Define a quadratic Poisson bracket by the formula 



Then the system can be written in the form Xi = {xi,H}, where H — 
The Louvillc intcgrability in the three-dimensional case can be easily established. 
In addition to the Hamiltonian function H, there exists a second intergral, in fact 
a Casimir F. The formula for this Casimir is given afterwards. We have to point 
out that in general algebraic intcgrability does not imply Liouville intcgrability and 
vice versa. 

In this paper we restrict our attention to the three dimensional case. For n = 3 
the system is defined by the matrix 



where a,b,c are constants. We use the notation {a,b,c) to denote this system. It 
turns out that the algebraically integrable Lotka-Volterra systems fall cither into 
two infinite families or four exceptional cases: 

Theorem 1. The Lotka-Volterra equations in three dimensions are algebraically 
complete integrable if and only if (a, 6, c) is in the class of 



[h) (1,0,1) 

(^3) (1,-1,1) 

(^4) (1,-1,2) 

(^6) (1,-2,3) 

Hx) (1,1, A) AeZ\0. 

ih) (i,h-m,m) MeR\0. 



We use the notation Zj to indicate that the system has an invariant of degree j 
or equivalently that the largest Kovalevski exponent is j. 



{■^z , -^j } — Clij XiXj , Z,J — 1,2,...,?!. 
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In Section 2 we give the basic definitions of weight-homogeneous vector fields 
and the Kowalevski matrix. The definition of Kowalevski exponents and relevant 
results follow the recent book [2] ■ See also the review article of Goriely [T^ where 
one can find many more properties of these exponents. In Section 3 we give some 
related properties of the Kowalevski exponents and a criterion of algebraic complete 
integrability. In Section 4 we define the three dimensional Lotka-Volterra systems 
and find necessary conditions for their algebraic integrability by analyzing the cor- 
responding Kowalevski exponents. In Section 5 we show that our classification is 
complete. Finally, in Section 6 we exclude any solutions that may exist due to 
higher order poles. 

2. Basic definitions 

We begin by defining what is a weight homogeneous polynomial. We follow the 
notation from [21 . 

Definition 1. A polynomial / G C [xi, X2, . . . , Xn] is called a weight-homogeneous 
polynomial of weight k with respect to a vector v = (ui, ^2, . . . , Vn) if 

f{t'"^Xi, . . .,f""Xn) = t^f {xi,X2, ■■■,Xn) ■ 

The vector v is called the weight vector. The Vi are all positive integers without a 
commont divisor. The weight k is denoted by w{f). 

Definition 2. A polynomial vector field on C", 

ii = fl{xi,X2, ■.■,x„) 

(2) ; 

Xji — fn {xi , X2 , ■ ■ ■ 7 ^n) 

is called a weight-homogeneous vector field of weight k (with respect to a weight 
vector v), if w(fi) — Vi + k = Tu{xi) + k for i = l,2,...,n. A weight-homogeneous 
vector field of weight 1 is called weight-homogeneous vector field. Furthermore, 
when all the weights are equal to 1, this is simply called a homogeneous vector 
field. 

Example 1. We consider the periodic 5-particle Kac-van Moerbeke lattice that is 
given by the quadratic vector field 

(3) Xi — Xi{xi—i — Xj^i), z = 1, . . . , 5, 

with Xi = Xi-f-<^ . This system has three independent constants of motion, 

Fi = Xl + X2 + X3 + X4+ X5, 

(4) F2 = X1X3 + X2X4 + X3X5 + X4X1 + X5X2, 

F3 = XiX2X3XiX5. 

Taking v = (1, 1, 1, 1, 1), becomes a homogeneous vector field and the weights 
of the integrals of motion are w{Fi)=l, w{F2) = 2 and w^F^) = 5. 

We now give the definition of an algebraically completely integrable system fol- 
lowing IT31[T3]. Note that this definition differs from the one given in [5]. However, 
the definition in [2] implies the definition of this paper. 

Definition 3. A vector field, 

il = fl{xi,X2, ...jXn) 



•^n — fn 1 5 -^2 5 ■ ■ ■ 5 ) i 
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is called an algebraically completely integrable system (a.c.i.) if its solution can 
he expressed as Laurent series 

oo 
fc=0 

where n ^ 1 of the coefficients are free parameters. 

2.1. Kowalevski Exponents. The following Proposition is important for two rea- 
sons. First, it gives an induction formula for finding the Laurent solution of a 
weight-homogeneous vector field and second it defines the Kowalevski exponents 
which is an important tool for our classification. 

Proposition 1. Suppose that we have a weight-homogeneous vector field on C" 
given by 

Xi — fi , . . . , Xji^ , i — 1,2,. ..,Ti, 

and suppose that 



(5) 



fc=0 



is a weight-homogeneous Laurent solution for this vector field. Then the leading 
coefficients, xf^\ satisfy the non linear algebraic equations 



(6) 



0, 
0, 



while the subsequent terms x^'^'^ satisfy 

(7) (fcId„-K;(a;(°)))x« 



wh 



(k) ^ 



ere x^ ' = 



( \ 

pends on the variables x^' 
matrix IC are given by 

(8) 



(0 



Xn^ with < I < k only. The elements of the n x n 



is a polynomial, which de- 



V Rn '' J 



dfj_ 

dxj 



ViSij , 



whe 



the Kronecker delta. 



Remark 1. The pole order vi of Xi in 0) coincides with the ith component of the 
weight vector. The number, Vi, is not necessarily the pole order of Xi because some 
of the x'^^^^ that can he calculated solving (0) may be equal to zero. 

Definition 4. 

The system 0) is called the indicial equation and its solution set is called the 
indicial locus and it is denoted by I. The n x n matrix JC, defined by 0), is called 
the Kowalevski matrix and its eigenvalues are called Kowalevski exponents (a 
terminology due to Yoshida). 



A necessary condition for algebraic integrability is that n — 1 eigenvalues of /C should 
be integers. It turns out that the last eigenvalue is always —1. The eigenvector 
that corresponds to —1 is also known. We have the following Proposition which 
can be found in [2]. 
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Proposition 2. 

For any m which belongs to the indicial locus X, except for the trivial element, the 
Kowalevski matrix /C(to) of a weight homogeneous vector field always has —1 as 
an eigenvalue. The corresponding eigenspace contains (wimi, . . . , w„m„)"^ as an 
eigenvector. 

3. Properties of Kowalevski exponents 

In this section we state some properties of Kowalevski exponents clarifying the 
connection with the degrees of the first integrals. We also give a necessary condition 
for a system to be algebraically completely integrable. The following results can be 
found in [H [IH [HI [HI El] . 

Theorem 2. // the weight-homogeneous system x = f{x) has k independent alge- 
braic first integrals Ii, . . . , Ik of weighted degrees di , . . . , dfe and Kowalevski expo- 
nents . . ■ , Pn, then there exists a fc x (n — 1) matrix N with integer entries, such 
that 

n 

'^Mij Pj = di, i^l,...,k. 

From this theorem wc have the two following corollaries: 
Corollary 1. 

// the Kowalevski exponents are 7i-independent, then there is no rational first inte- 
grals. 

Corollary 2. // the Kowalevski exponents are ISI -independent, then there is no 
polynomial first integrals. 

We also have the following theorem, see [T^ [T51 [^ . 

Theorem 3. Suppose that the system (QJ) possesses a homogeneous first integral 
Fm of degree m. Then there exists a set of non-negative integers fe, • • ■ , fc„ such 
that 

n 

fcj Pj = m A:2 + /ca + • • • + fc„ < m . 

The next theorem which can be found in 2] gives us a necessary condition for a 
system to be algebraically integrable. This criterion can be checked easily simply 
by computing the Kowalevski exponents. 

Theorem 4. 

Let pi ~ —1. A necessary condition for a system of the form ^ to be algebraically 
completely integrable is that all the Kowalevski exponents p2, . . . , Pn should be inte- 
gers. 

4. LOTKA-VOLTERRA SYSTEMS 

4.1. Hamiltonian formulation. Consider a Lotka-Volterra system of the form 

n 

(9) Xj = ^ ttjkXjXk, for j = 1, 2, . . . , n, 

fe=i 

where the matrix A = (aij) is constant and skew symmetric. 
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There is a sympletic realization of the system which goes back to Volterra. In 
other words a projection from R^" R" from a symplectic space to a Poisson 
space. Volterra defined the variables 

qi{t) = / Ui{s)ds 

JQ 

and 

1 " 

p.i{t) = ln((7j) - 2 X! 

k=l 

for i = 1,2, ... ,n. Now the number of variables is doubled and Volterra's transfor- 
mation is given explicitly by 

X, = eP'+^ ^''-'=1 for i = l,2,...,n. 

The Hamiltonian in these coordinates becomes 

n n n 

H = j2^^ = J2'i' = J2 ^"'^^ ■ 

2 — 1 2—1 2—1 

The equations © can be written in Hamiltonian form 

If ={*'^}' 

-If ={P^,H}, 

i = 1,2, . . . ,n, and the bracket {•, •} is the standard symplectic bracket on R^": 

r 1 r f 1 1 if * = j ■ • 1 f> 
{qt,Pj} = S^J = i , t,.] = l,2,...,n ; 

all other brackets are zero. The corresponding Poisson bracket in x coordinates is 
quadratic 

{ ; -^j } — O^ij -^j 1 ^: J — 1 , 2, . . . , TL. 

Equations © in a; coordinates are obtained by using this Poisson bracket and the 
Hamiltonian, H = xi + X2 + ■ ■ ■ + Xn- 

4.2. The three-dimensional case. 

In this paper we restrict our attention to the three dimensional case. For n = 3 
the system is defined by the matrix 

/ a b 

(10) A = [ -a c 

\ -b -c 

where a, b, c are real constants. 

Using equations ([8]) we obtain the Kowalevski matrix 

ax^^ + bx^^ + 1 ax^^ bxf*^ 

(11) I -axf _a4°)+cxf + l cxf 
—bx^^^ —cx^^^ —bxf'^ — cx^2^ + 1 

where x^^^ = (^xf'\ x^2 \ x's^^ is an element of the indicial locus, i.e., a 
solution of the simultaneous equation which in this case is written as 

(12) xf^-axf^xf^ 
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In Table 1 we list the corresponding Kowalevski exponents for each element 
of the indicial locus. 



Vector x^^^ 


Kowalevski 
exponents 


Vector x^'^) 


Kowalevski 
exponents 


(0,0,0) 


1,1,1 

1 1 a—b+c 


(^ri,0) 


_ll a-b+c 
' ' c 

_ll a-b+c 
' ' a 



Table 1 . Kowalevski exponents of 3x3 Lotka-Volterra equations 



A necessary condition of algebraic integrability is that all the Kowalevski ex- 
ponents must be integers. So we have to solve the simultaneous Diophantine 
equations 

a-b + c , a-b + c , a-b + c 

(13) = fci, = k2, 7 = fes, 

a c b 

where ki, k2, k^ € Z. The case b = a + c for which ki = k2 = k^ = is 

investigated below. Solving (fT3]) we find that 

t^i4j K3 — ) ^ k2^ L _ fci+fc2-fclfc2 



hk2 -ki -k2 1 « = 17^, 



k 



c 



(15) 



(16) ki 



hks jfe = -|«> c=^i±^|^a 

k,ks-k,-ks' \ a = -^^b, c= 

k2k3 /^ = -!tc, a = ^^^^^f^c 

k2ks-k2-k3' ] c=-H, a=^2±^|^6 



fc2 ~ fc2 

We assume first, that the Kowalevski exponents are not zero. We examine 
the solution 

kik2 , ki + k2-kik2 ki 

fcifc2 - ki - k2 k2 k2 

We determine the values of ki and A;2 so that the fraction, 

, kik2 

kik2 - ki - k2 

is an integer. We first consider the case kik2 — /ci — ^2 7^ 0. 

Case I Assume positive values for both ki and /c2- 
Since 

kik2 = 1 + ki+k2 
kik2—k\—k2 k\k2—ki—k2 

it is enough to to solve the Diophantine equation 

x + y 

= z 

xy — X — y 

for X, y positive integers and z G Z. 
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Lemma 1. Let x,y & Z+ with x <y. Then 

x + y ^2 



xy — X — y 



if and only if {x,y) is one of the following: (1,A), A G Z"*", (2,3), (2,4), 
(2,6), (3,3), (3,6), (4,4). 

Proof. Since 

xy — X — y<x + y 

we have 

xy < 2{x + y) <4y . 

Since y we get x < 4. Therefore x = 1,2,3,4. We examine each case 
separately. 

• If a; = 1 

x+y 1 +v , „ 

^ = = -l-yeZ. 

xy — X — y —1 

Therefore (1, A), A G Z+ is always a solution. 

• Suppose X = 2. Then 

x+y 2+y ^ , 4 



xy — X — y y — 2 y — 2 

should be an integer. Therefore y — 2 = ±1,±2,±4. We obtain the 
solutions (2,3), (2,4) and (2,6). 

• Suppose X = 3. Then 

x + y ^ y + 3 
xy — X — y 2y — 3 

should be an integer. Therefore 

2y - 3 < y + 3 

and we obtain y < 6. We obtain the solutions (3, 3) and (3, 6). 

• Suppose X = A. Then 

x + y y + 4 



xy — X — y 4y — 4 
should be an integer. Therefore 

3y - 4 < y + 4 
and we obtain y < 4. We obtain the solution (4,4). 

Of course, since the fraction 

x + y 



□ 



xy — a; — y 

is symmetric with respect to x and y, we easily obtain all solutions in positive 
integers. 

We summarize: 
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(18) For l<ki<k2 < 



' ki = i =^ ^2 = A e z+ 

ki=2 =^ /C2 e {3, 4, 6} 

ki = 3 => /C2e{3,6} 

ki = i =^ k2 = 4. 



Note that the case fci = 3, A;2 = 3 implies k^ = 3 and we obtain the periodic 
KM-system (1,-1,1). 

Case II 

Suppose one of them, say ki, is positive while the other, k2, is negative. Let 
/c2 = —X, X > 0. Then 

—kix kix , X — ki 

ks = -r-— = 1 = 1 + 



-kix — ki + X kix + ki — X kix + ki — x 

It is enough to to solve the Diophantine equation 

x-y 



= z 



xy + y - X 
for x, y positive integers and 2; G Z. 

Lemma 2. Let x,y E Z+. Then 

x-y 



€ Z 



xy + y - x 

if and only if (x, y) is of the form (A, 1) or (A, A) with A G Z"*". 
Proof. li y = 1 then 

X — V „ 



xy + y — X 

Therefore a pair of the form (A, 1) is always a solution. 
Assume y > 1. We note that 

xy x-y 



xy + y — X xy + y — X 

and therefore xy + y — x < xy implies y < x. li x = y then our fraction is 
clearly an integer. On the other hand, if y < x, then the fraction 

X — y . „ 

— # Z 

xy-\-y — X 

since 

(y-l)x + y>x + y>a;-y. 

□ 

If fci = 1 then k2 = —A and k^, = A. Similarly, \ik\ = \ then k2 = —A and 
^3 = 1. The two cases are isomorphic and correspond to Case 5 in Table 3. 

Case III 

If we take negative values for both ki and k2, then 

, xy x + y 

ks = = 1 , 

xy + x + y xy + x + y 
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where ki = —x and ^2 = —y with x,y > 0. We have that xy > and 
xy + X + y > so that the Kowalevski exponent is an integer if 

xy + x + y<x + y 

which imphes xy < 0, a. contradiction. Therefore, in this case cannot be 
an integer. 

This completes the analysis of the case kik2 — ki — k2 ^ 0- 
Now suppose kik2 — fci — ^2 = 0. 

In this case we have ki + k2 = kik2 and obviously (since we assume non- 
zero Kowalevski exponents) we must have ki = k2 = 2. We easily obtain 
a = c and 5 = 0. This system is equivalent to the open KM-system (also 
known as the Volterra lattice). This is Case 1 in Table 3. 

This concludes our analysis. We have obtained necessary conditions for 
the algebraic integrability of Lotka- Volterra systems in three dimensions and 
the results are summarized in Table 2. In Table 3 we also include the case 
of a zero exponent i.e. b = a + c. Note that the case b = a + c which is 
equivalent to (1, 1 + /i,/i) for fi e R \ was also considered in [7] from a 
different point of view. 

4.3. Equivalence. 

In order to have a more compact classification, we define an equivalence 
between two Lotka- Volterra systems. To begin with, common factors can 
be removed. In other words, suppose that matrix A = (uij) in ([9]) has a 
common factor a. Precisely, if 

aij = Cija, where Cij G R, i, j = 1,2, . . . ,n, 

then the Lotka- Volterra system ([9]) can be simplified to 

n 

U{ = ^ ^ CijUiUj , Z = 1,2..., Tl, 

i=i 

using the transformation 

Ui = axi, i = 1,2 ... ,n. 

More generally, we consider two systems to be isomorphic if there ex- 
ists an invertible linear transformation mapping one to the other. Special 
cases of isomorphic systems are those that are obtained from a given system 
by applying a permutation of the coordinates. Let a G Sn, and define a 
transformation 

— >x^(^i), i = l,2,...,n. 

The transformed system is then considered equivalent to the original sys- 
tem. We illustrate with an example for n = 3. 

Example 2. We prove that the system 



Xl = axiX2 — f 3^1X3 Xl = 3xiX2 — X1X3 

±2 = -axiX2 + ^^2X3 Ui = axi ' ±2 = -3x1X2 + 2x2X3 

X3 = |X1X3 - ^X2X3 X3 = X1X3 - 2X2X3 

is isomorphic to the system 
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xi = axiX2 — 2axiX3 xi = xiX2 — 2x1X3 

X2 = —axiX2 + 2>ax2X'i Ui = axi ' X2 = —X1X2 + 8x2X3 
= 20x1X3 — 80X2X3 X3 = 2x1X3 — 8X2X3 

Applying a = (1 8 2) we have that 

= =xz= X1X3 - 2x2x3 = X2X1 - 2X3X1 

X2 = X^(2) =Xi= 8X1X2 - X1X3 = 8X2X3 - X2X1 , 

^3 = 2:^(3) = X2= -8x1X2 + 2x2X3 = -8X2X3 + 2X3X1 
which is the second vector field. 
Example 3. Note that the system 

Xl = -X2X3 

X2 = 2:2X3 

X3 = X1X3 - X2X3 - x| 

is equivalent to the open KM-system (1,0,1) under the transformation 

(Xi, X2, X3) (X2 + X3, Xl, X2) 

hut it is not a Lotka-Volterra system. 

In Table 2 we display the different values of (a, b, c) of the solutions ([H]), (fT5]) 
and (jl6p of the simultaneous equations (jl8p which ensure integer Kowalevski 
exponents for the Lotka-Volterra system in three dimensions. We also list 
the elements of the symmetric group S3 which realize the isomorphism. Note 
that A € Z \ 0. The final six non- isomorphic systems are displayed in Table 
8. 



Vector (a, b, 


c) 


Kowalevski 
exponents 


a 


(a, a, Xa) 
(a, Xa, —a) 


-1,1,1 
-1,1,A 
-1,1, -A 


0- = (1 3) 
0- = (2 3) 


(a, —a, 2a) 
(a, —2a, a) 


-1,1,2 
-1,1,4 


0- = (1 3) 
0- = (1 3 2) 


(0, —a, a) 


-1,1,8 




( 
( 


^ a 

a, 

k"' 2 ' 2; 

^ a 

k"' 2 ' 2 ^ 

a, —2a, 8a 
a, —3a, 2a 


) 
) 


-1,1,2 
-1,1,3 
-1,1,6 


0- = (1 3 2) 
0- = (1 3) 
0- = (1 2 3) 
0- = (2 3) 

0- = (1 2) 


(a, 0, a) 
(a, -a,0) 
(0,6,-6) 


-1,1,2 


0- = (1 3 2) 
0- = (1 2 3) 



Table 2. Systems with integer Kowalevski exponents 
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Example 4. The periodic KM system ([15]) in three dimensions is the sys- 
tem 



^ ^ CLjjXjXj J i — 1,2,3, 



Xi 

i=l 

where A is the 3x3 skew- symmetric matrix 



A 




This system is a special case of the system (0j where (a, b, c) = (—1, 1, — !)• 
This is Case 2 in Table 3. The Kowalevski exponents of this system are 
— 1, 1,3. The system can be written in the Lax-pair form L = [L,B], where 

/ xi 1 \ / 

L = i 1 X2 , B = I X2X3 

\ X3 1 y \ X1X3 

We have the constants of motion 




Hk = trace l^L" j , k=l,2,... 
The functions 

H2 = Xi + X2 + X3 
H3 = 1 + X1X2X3 

are independent constants of motion in involution with respect to the Poisson 
bracket 

(0 -X1X2 X1X3 
X1X2 -X2X3 
-X1X3 X2X3 

We note that the positive Kowalevski exponents, 1 and 3, correspond to the 
degrees of the constants of motion. 

We have to point out that ah Lotka-Volterra systems in three dimensions are 
integrable in the sense of Liouville since there exist two constants of motion 
which are independent and in involution. The function 

-f/' = Xi + X2 + X3 

is the Hamiltonian for these systems using the quadratic Poisson bracket 

(0 axiX2 6x1X3 

— axiX2 CX2X3 
—6x1X3 — CX2X3 

The equations of motion can be written in Hamiltonian form 

— {xj, i?}, z = 1, 2, 3 . 

The second constant of motion, independent of H always exists. It is 
straightforward to check that the function 

F = xlx^'^x^ 

is always a Casimir. Therefore the system is Liouville integrable for any 
value of a, b, c. This is not the case if n > 4. 
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5. Free Parameters 

We would like to classify the algebraic completely integrable Lotka-Volterra 
equations in three dimensions. In order to use Proposition [T] we have to as- 
sume Laurent solutions of the form 

(19) ;r.(^) = _^:cf)^^ i = l,2,...,n, 

k=0 

where Vi are the components of the weight vector v that makes the vector 
field 

Xi — fi{x\i . . . 1 Xn)i i — 1)2,... 

to be weight homogeneous. In our case the weight vector is v = (1, 1, 1). 

To make sure that our classification is complete we must check that each 
system obtained by imposing these necessary conditions is indeed a.c.i. This 
means that the Laurent series of the solutions xi, X2 and X3 must have 
n — 1 = 2 free parameters. Using the results in [2j, the free parameters 
appear in a finite number of steps of calculation. The first thing to do is to 
substitute (fT9]) into equations ([9]). After that we equate the coefficients of 
t^. We have already equated the coefficients of t~^^~^ by solving the indicial 
equation to find xf^\ Then we call Step m {m € N) when we equate the 
coefficients of ^-^^-1+™- to find x["^\ According to [2] all the free parameters 
appear in the first kp Steps, where kp is the largest (positive) Kowalevski 
exponent of the system. The calculations are straightforward and we omit 
the details. 



Vector 


Kowalevski 


Free parameters 


Degree of 




exponents 






invariant 


(0,0, a) 


-1,1,2 


X3 , 




2 


(a, —a, a) 


-1,1,3 


■^3 ' 


^(3) 
X3 


3 


(a, — 2 , 2) 


-1,1,2 
-1,1,4 


X3 , 


(2) 


4 




-1,1,2 








(a, —2a, 3a) 


-1,1,3 
-1,1,6 


X3 , 


X3 


6 




-1,1,1 










-1,1, A 
-1,1, -A 


■^1 , 


X2 


A 


(a, a + c, c) 


-1,1,0 


x^ , 


AO) 






Table 3. Free parameters for the algebraic complete integrability 



of each system 



All the systems that we have obtained turn out to be a.c.i. We summarize 
the results in Table 3 where we display the 2 free parameters in each case. 
We note that the six Cases of Theorem 1 are non-isomorphic by examining 
the degree of the Casimir. 
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6. Higher Order Poles 

In our classification, using the Kowalevski exponents, we assume that the 
order of the poles agrees with the components of the weight vector, in our 
case all equal to 1. We have to exclude the possibility of missing some cases 
due to solutions with higher order poles. We show that no such new cases 
appear. 

Suppose that the Laurent solution of the system is 



^iW = FrEr=o^i with4"V0, 
(20) X2{t) = ^ Er=o with / 0, 

^3{t) = ^j:r=o^3 ^t^ withxf/o. 

If 1^1 , z/2 , z/3 < 1, then these systems have been already investigated using 
Proposition [H On the other hand, keeping in mind that H = xi + X2 + X3 
is always a constant of motion, we end-up with the following four cases to 
consider: 

(i) z^i = f 2 = > 1 and ^3 < u, or 

(ii) vi = 1^3 = I' > 1 and 1^2 < or 

(iii) 1^2 = 1^3 = > i and z^i < u, or 

(iv) Vi = 1^2 = = > ^■ 

Recall that equations ([9]) in three dimensions are: 



(21) 
(22) 
(23) 



Xl 
X2 



±3 = -6x1X3 

We examine each of the four cases: 



axiX2 + 6x1X3, 
axiX2 + CX2X3, 
CX2X3. 



(0 



1^1 



^2 = f > 1 and 1/3 < p 



Since ui = U2 = v and using the fact that H 
constant of motion we have that 



Xl + X2 + X3 is a 



a / . 



We also note that u + V3 < 2v and x^ X2 7^ 0. Equating the 
coefficients of t^"" of the LHS and RHS of ^ or ([22]), we are led to 
a xf^x^^ = 0. Therefore a = 0. 

As we know that U3 + 1 < 1^ + 1/3, the coefficient of t'^+^'s Qf ^j-^g 
RHS of ([231) must be equal to zero. So 

but Xg'^^ / and x'^'^ = —x^^^ ^ 0; therefore b = c. 
If 6 = 0, then from ^ and ^ we have that 

Xl = ±2 = =^ Xl, X2 are constant functions, 

that is a contradiction because z^i = z^2 = > 1- 
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(24) 
(25) 



If b and c are non-zero, then the equations (j21|) and (j22p become 

Xl = 6X1X3, 
X2 = 6X2X3. 

Using equation (|2T]) we obtain 
1/ + 1 = 1/ + 

therefore 2/3 = 1, since x^'^^Xg'^^ / and x^2^x^^^ 7^ 0. 
It follows from jMl) and (l25|) that 

Xl X2 



■ = 6x3 = 

Xl X2 

However, we know that 



Xl = KX2, K is a constant. 



.(0) 



„(0) 
"2 



Equation (p3|) becomes 
X3 = -6(-X2)x3 - 6X2X3 



Xl 



-X2. 







This is a contradiction since 1^3 = 1 and xr^ / 0. 



X3 = c, c is a constant. 

(0) 



(u) 
{in) 
(iv) 



z/3 = z/ > 1 and V2 < ^ 



It leads to a contradiction, as in case {i). 



1^2 = = ^ > ^ and Ui < p 



It leads to a contradiction, as in case (i) 



/^i = Z/2 



V > 1 



In this case, for i = 1, 2, 3, 



Xi 



A;=0 



we have that the degrees of the leading term of the LHS of the equa- 
tions (pTI) . (f22|) and (i23]l are equal to u + 1, but the degrees of the 
leading term RHS of these equations are equal to 2i/ and so the 
coefficients of of the RHS of these equations must be zero for 
= 2, 3, . . . , z^. 

The coefficients of k = 1,2, . . . ,1^, are given by the sums 



v—k 



(26) 



Si,k = ^ xf'^'^f^k^ foi^ i = 1, 2, 3, 



A=0 



where 



u 



u. 



(A) 

1, fc 

(A) 

2, fc 
(A) 

3, fc 



ax 



(u-k-x) .{y-k-x) 

2 ^ "-^3 ' 



-ax 



-bx 



(u-k-X) 
1 

(u-k-X) 



+ CX 



CX 



(v-k-X) 
3 ' 
(i/-fc-A) 
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Note that 



= ^Ij "* k + X = j + m. 



In addition 

Si,k = 0, for i = 1, 2, 3 and A; = 2, 3, . . . , i/. 
For k = n sum (1261) becomes 



smce 7^0. 

For k = v — 1 we have that 

.(o),,(o) , ™a)„,{o) _ (0) (0) _„ 

=^ = ^ because x^- 7^ U. 

Let m G {1, 2, . . . , 1/ — 1} and assume that u[^l = for /c > m. 
For k = m we have that 

u—m u—m 

„W„,(A) _ ^(0)^,(0) , ^(A) (A) 



c. _ V x^^K^^^ - X^^K^^^ + V 
Ji,m — Xj u,^ ,^ — ttj „^ -r 



A=0 A=l 

u—m 

i,m+A -^i "'i,m- 



A=l 

Since Si^m = for m > 1 and, since xf"* 7^ 0, then ufl^ = 0. 
Now we equate the coefficients of ^^qrr on both sides of the equa- 
tions (IlID-dSiD to obtain 

C (0) (0) (0) 

Therefore, v + ufl = 0. 
Therefore we have that 

— axj 
-hx[ 

These simultaneous equations have solutions only if 

h = a + c. 

If a = 0, Then b = c (obviously 6 = c 7^ 0). Then the system is 
isomorphic to the following (0, 1, 1) system: 

Xl = X1X3 
X2 = X2X3 

X3 = -X1X3 - X2X3 

Equating the coefficients of t^^'^ {v > 1) in the first and second 
equations we have that 

™(0)^(0) _ ^(0) (0) _ ^ 



2 


+ 




= 




1 


+ 


(u- 

cx\ 


= 




■2 




{u- 
CX3 


= 


— V. 
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This is impossible because x-^-* / 0, for i = 1,2, 3. 

The same happens if be = 0. So in the following calculations we 
assume that abc ^ 0. 

We will show that there exists no such solution with v >2. Since 
b = a + c and the function H = xi + X2 + is a constant of motion, 
the Lotka-Volterra equations in three dimensions can be written in 
the form 

xi = akxi — ax\ + cxix^, 
(29) ±2 = -xi-is, 

xs = —ckx3 + cx| — axixs, 

where k is the constant value of the function H. It is straightforward 
to see that if A: ^ 0, then the solution is 



kCie^''^ kae-^^^ 

^1 = 77 



(^0) X2 = k - xi - X3 = - , ^^'Lcfct — ?r- 

Obviously C2 / 0. The pole t* satisfies 

Cie"^^* + ae'"^^* - C2 = ^ C2 = Cie"^^* + ae'"^^* + 
Hence using De 1' Hopital Rule we are led to the fact that 

C2 

lirn it - U)x2(t) = , , — r^. 

Since the pole order is greater than 1, we have that 
lim [t — t*)x2(i) = 00. 

t — ^■t * 

Therefore 

The solution ([HU]) possesses only one arbitrary constant k, but we 
need 2. 

Now if /c = the solutions of (1^^ are 



X3(t)=0, xi{t) ^ 



at + Ci 
or 

, , Ci — c , . 1 



a{Cit + C2)' ' Cit + C2 ■ 
Both solutions lead to a contradiction since the pole order of xi 
and X3 is assumed to be greater than 1. 

Therefore the case z^i = 1^2 = = > 1 does not give us any new 
algebraically integrable systems. The conclusion is that the case 
6 = a + c is algebraically integrable only when z/i = ^2 = fs = 1. 
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